Deformation and its influence on K isomerism in neutron-rich Hf nuclei 
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We investigate the influence of deformation on the possible occurrence of long-lived K isomers in 
Hf isotopes around TV = 116, using configuration-constrained calculations of potential-energy sur- 
faces. Despite having reduced shape elongation, the multi-quasiparticle states in i^^.issjjj- j-gu-^^in 
moderately robust against triaxial distortion, supporting the long expected occurrence of exception- 
ally long-lived isomers. The calculations are compared with available experimental data. 

PACS numbers: 21.10.-k, 21.60.-n, 23.20.Lv, 23.35.+g 
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Atomic nuclei have special stability at some proton or 
neutron numbers named magic numbers. Analogously, 
there exist some K-vaagic (K is the total angular momen- 
tum projection onto the symmetry axis of a deformed nu- 
cleus) numbers associated with exceptionally long-lived 
high- ivT isomers P, Q • Such isomers are formed through 
broken-pair or multi-quasiparticle (multi-qp) excitations. 
The long half-life usually originates from a combination 
of relatively-low excitation energy, high K value and well- 
deformed axially-symmetric shape. These factors favor 
the conservation of the K quantum number and hence 
the hindrance of decay to low- if states. The best known 
-ftT-magic number is probably Z = 72 where K'^ = 8~ iso- 
mers with the configuration 7r2{7/2+[404] (g) 9/2-[514]} 
systematically occur in ^''^"i^^Hf, with half-lives ranging 
from nanoseconds to hours 0, The neutron num- 
ber N — 106 is another well-known ii'-magic number. 
Here > 1-/^8 = 8~ isomers with the configura- 
tion z/2{7/2-[514] (g) 9/2+[624]} are observed from ^^''Er 
to i^^Pb 1^. These scarce long chains of isomers provide 
unique opportunities for systematic study. The doubly 
if-magic nucleus ^^^Hf manifests itself with the occur- 
rence of a 31-yr K"^ = 16"'" isomer based on the coupling 
I'g- <8) TTg. . The possible stimulated release of the stored 
energy (2.4 MeV) in the isomer has been associated with 
proposals for a clean reservoir of nuclear energy and 7-ray 
lasers ^MM- 

Exceptionally long-lived isomers have been predicted 
to also occur in neutron-rich Hf isotopes due to favorable 
conditions, i.e., low energy and high K value 0, 0, 
With a = IS"*" isomeric state from the configuration 
:/2g„{9/2-[505](8)ll/2"^[615]}®7r|_, ^^SHf seems to be a 
doubly if-magic nucleus [2!]. Experimentally, projectile- 
fragmentation ^] and deep-inelastic [IQ] reactions have 
been employed to explore this region, such as the obser- 
vation of a 240-/iS isomer with the configuration in 
major advance was achieved recently us- 
ing storage-ring mass measurements f4| , where long-lived 
isomers were observed in neutron-rich Hf isotopes up to 
i^^Hf, with iV = 114. Nevertheless, the predicted doubly 
JT-magic character of ^^^Hf, with N = 116, remains to 
be confirmed. 



However, another important factor, deformation, for 
the emergence of K isomers could be much different be- 
tween i^^Hf and the nuclei around ^^^Hf. With neu- 
tron number N = 106 at the mid-shell between spherical 
magic numbers 82 and 126, ^^^Hf is well deformed, having 
a shape that favors the formation of long-lived isomers. 
In contrast, a reduced shape elongation is expected in 
i^^Hf because the neutron number N — 116 is far away 
from the mid-shell. Generally, a smaller /32 deformation 
is more susceptible to shape fluctuation towards triaxial- 
ity that induces K mixing. It is still an open question to 
what extent the deformation influences the occurrence of 
exceptionally long-lived A'-isomers in neutron-rich Hf iso- 
topes. In the present work, we investigate these high- if 
isomeric states using configuration-constrained calcula- 
tions of potential-energy surfaces (PES) Jd.,J^], a model 
that properly treats the deformations of multi-qp states. 

In our model, the single-particle levels are obtained 
from the deformed Woods-Saxon potential with the uni- 
versal parameter set [13, To avoid the spurious 
pairing phase transition encountered in the usual BCS 
approach, we use the Lipkin-Nogami (LN) treatment of 
pairing [l5j, with pairing strength G determined by the 
average-gap method [l3|- The configuration energy [ll| 
in the LN approach can be written as 
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where S is the proton or neutron seniority for a spec- 
ified configuration (i.e., the number of blocked orbits 
with index kj), and N is the proton or neutron num- 
ber. The orbits with index kj are blocked at each point 
of the (/32,7,^4) [H or (h^PA^Pa) [12^ deformation lat- 
tice, which is the so-called configuration constraint. This 
is achieved by calculating and identifying the average 
Nilsson quantum numbers for every orbit involved in 
a configuration. The total energy of a configuration 
consists of a macroscopic part which is obtained from 
the standard liquid-drop model [l7j and a microscopic 
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TABLE I: Neutron configurations of tiie i/^(g)7r^_ {7/2+ [404] 
9/2~[514]} isomeric states. 
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FIG. 1: Excitation energies (a), hexadecapole deformations 
(b), and quadrupole deformations (c) for multi-qp isomeric 
states in Hf isotopes. Open and filled swnbols show the cal- 
culated and experimental values [1, [EISli respectively. The 
data for the 2-qp state in ^^*Hf (A'^ = 106) corresponds to an 
unperturbed 7rg_ state (see text). Note that no experimental 
configuration is available for the nuclei with N > 112. The 
deformations of 4-qp states (not shown for clarity) are similar 
to those of the corresponding 2-qp states. For all the multi-qp 
states, 7 « and Pa ~ 0. The same configurations are con- 
nected by lines to guide the eye. See Table |I] for the neutron 
configurations of the 4-qp isomeric states. 



part resulting from the Strutinsky shell correction |18| . 
SEshcii — -Eln — Egtmt- The minimum of the PES gives 
the energy, deformation and pairing property of a multi- 
qp state. 

Figure [1] displays the calculated excitation energies 
and deformations of high-K isomeric states in even-A 
i76-i92jj£ isotopes. Our calculations satisfactorily repro- 
duce the excitation energies of the tt^ 8~ isomers observed 
in 176-182 jj£ ^ should be noted that the observed 4-s 
isomer in ^^^Hf has a mixed configuration of 64% i^g_ and 
36% TTg. [13] , resulting from the interaction between the 
two close-lying and tt? configurations. The exper- 
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imental excitation energy of the unperturbed 7r|_ state, 
1.35 MeV [5], is well reproduced by our calculation of 1.35 
MeV. The calculation of the 8" state in ^^^^Hf gives an 
excitation energy of 1.23 MeV that is in good agreement 
with the measured value of 1.26 MeV [4i]. This indicates 
the association of the isomer with the tt^_ configuration, 
which is consistent with the suggestion in Ref . [J] . 

The calculated excitation energy of the K"^ = 8" state 
suddenly increases at iV = 118, as shown in Fig. [TJ Here 
there occurs a ground-state phase transition from pro- 
late shape to oblate shape. Our prediction is consistent 
with other calculations (see, e.g., Refs. [Ill, [l^])- In con- 
trast, MoUer et al. predict that ^^°Hf is still prolate de- 
formed [11] , indicating that the prediction of the neutron 
number for the prolate-to-oblate transition is model de- 
pendent. The 8" state in ^^°Hf is calculated (see Fig.[T]) 
to be at an excitation energy of 1.86 MeV, but has a pro- 
late shape, distinct from the oblate shape of the ground 
state. Also the calculated (34 deformations are much dif- 
ferent between the isomeric state and the ground state. 
A combination of K isomerism and shape isomerism may 
happen in ^^"Hf, analogous to that predicted for the i/^ 
8" isomer in ^^^Hg |24|. Whether or not the higher exci- 
tation energies and smaller j32 deformations predicted for 
the more neutron-rich Hf isotopes will lead to the termi- 
nation of the K"^ — 8" isomer chain at A^ sa 118 remains 
to be determined. 

Figure [T] also shows the lowest-lying 4-qp high-if 
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Hf isotopes that are formed 
through the coupling of i/"^ (see Table |T| with 7rg_ . For 
^**Hf, the K'^ = 18+ state has calculated excitation en- 
ergy close to that of the = 16" state. The calculated 
excitation energies are reasonably in accord with experi- 
mental data of ™-^^^m. The observed 2.48-MeV isomer 
in ^*^Hf is likely the lowest-lying 4-qp 15+ state. The cal- 
culated excitation energy of the 17+ state is closer to the 
experimental data, but it can barely be long-lived be- 
cause of the possibility to decay to the 15+ state. The 
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FIG. 2: Excitation energy relative to a rigid rotor for the 
lowest-lying 4-qp isomeric states in even-j4 Hf isotopes. Cal- 
culated and experimental values [1, 13] are represented by 
lines and squares, respectively. The data point for ^*®Hf is 
not shown as the K value is unclear. Note that the 16~ state 
in ^**Hf has larger Ex — Er value than the 18^ state because 
of the smaller angular momentum. 



newly observed isomer in -'^^^Hf was tentatively associated 
with the 17+ configuration in Ref. ^J] based on blocked- 
Nilsson-|-LN calculations with fixed deformations. That 
significantly underestimated the measured excitation en- 
ergy. In our calculations, the inclusion of deformation 
effects does not remove this discrepancy. 

The discrepancy in ^^^Hf could lie in both experiment 
and theory. Experimentally, the data for the 4-qp isomer 
in ^^^Hf, obtained from just two ions, need to be con- 
firmed. Furthermore, the mass measurements do not give 
a complete picture of the low-lying intrinsic states and do 
not provide spectroscopic information. Theoretically, our 
calculations do not include the residual spin-spin interac- 
tion [7| that can change the calculated excitation energy 
by up to several hundred keV. (We note that all the neu- 
tron configurations listed in Table U have energetically fa- 
vored spin-antiparallel couplings except for the 17+ state 
in ^^^Hf.) In addition, the Woods-Saxon single-particle 
levels could become less accurate when moving away from 
the stability line. It has been remarked previously that 
the observed 1-qp states in odd-A W nuclei approaching 

= 114 begin to deviate considerably from the predic- 
tions of the Nilsson model [2^ . Our model cannot gen- 
erate the correct state orderings in ^^'^W either, which 
implies a decreasing validity of the present Woods-Saxon 
potential parameters in the neutron-rich nuclei. These 
parameters are determined by fitting the properties of 
nuclei near the stability line. Finally, the pairing strength 
could be another factor influencing the calculated exci- 
tation energy. An enhancement of about 10% (relative 
to the value obtained from the average-gap method) for 
both protons and neutrons is obtained for i^°W by fitting 
experimental odd-even mass differences (26j . This incre- 
ment in pairing strength can result in an increase of the 
calculated excitation energy by « 400 keV for the 
state. Unfortunately, such a fitting process is not avail- 




-20 -10 10 20 -20 -10 10 20 



Y(deg) Y(deg) 

FIG. 3: (Color online) Calculated potential-energy curves as 
a function of 7 deformation. The energy is minimized with 
respect to /32 and fii deformations. The true minimum of the 
oblate-deformed ground state in ^^"Hf is at 7 = —60° in the 
Lund convention. 

able for the neutron-rich Hf nuclei due to the absence of 
mass data required to calculate the odd-even mass dif- 
ferences. The above-mentioned problems in residual in- 
teraction, single-particle energy and pairing strength still 
await to be properly treated to obtain a more complete 
and accurate description of the neutron-rich Hf nuclei. 

In Fig. [5J we plot the excitation energy relative to a 
rigid rotor for the 4-qp isomeric states. It has been found 
that the degree of K hindrance increases in a simple man- 
ner with decreasing excitation energy relative to a rigid 
rotor [H [3, This is because a relatively- low exci- 

tation energy is often associated with low statistical K 
mixing due to the low level density [27j , favoring the con- 
servation of the K quantum number. Figure [2] shows two 
minima at iV = 106 and N « 116, which is mainly as- 
cribed to the closeness of the two single-particle orbits 
occupied by the unpaired neutrons and coupled to high 
K. The minimum at w 116 is even lower than that at 
the X-magic number N — 106 where the well-known 31- 
yr K"^ = 16+ isomer occurs. This, although consistent 
with the prediction in Ref. 0, needs to be confirmed by 
further work because of the uncertainty in our calcula- 
tions of excitation energy around = 116 (see above). 
It is the remarkably low energy together with the high 
K value that indicates the possible emergence of excep- 
tionally long-lived isomers [1, 0, 0] ■ 

Nevertheless, these isomeric states have smaller P2 de- 
formations than ^^^Hf. It is usually easier for a less elon- 
gated shape to fluctuate towards 7 distortion. Axially 
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asymmetric deformation is intimately involved with K 
mixing, because there is then no axis of symmetry on 
which to define a K projection. Some high-_?C isomers 
with very low hindrance in decay have been found to 
have 7 deformations (see, e.g., Ref. [HI). Even though 
a K isomer has axially symmetric shape, it can decay 
by tunnelling through triaxial shapes [28|. The softness 
against 7 deformation facilitates such a decay mode. To 
see the degree of 7 softness accompanied by small /32 
deformation, we calculate the potential-energy curves as 
a function of 7 deformation for Hf isotopes, as shown 
in Fig. [31 The doubly if-magic nucleus ^^*Hf can be 
treated as a benchmark, to which the neutron-rich nu- 
clei are compared. The qualitative comparison can be 
barely influenced by the potential inaccuracy in pairing 
strengths and single-particle levels. It can be seen that 
the ground states of neutron-rich Hf isotopes become in- 
creasingly soft against triaxial distortion with decreas- 
ing 132 deformations. With \P2\ < 0.2, iss^iaoHf have 
remarkably flat ground-state potential wells. However, 
the multi-qp states in i86,i88jjf show enhanced stiffness in 
the potential wells compared to the corresponding ground 
states. The isomeric states in ^**Hf have potential wells 
not much different from those in ^*^Hf where a long-lived 
isomer has been observed [3] , though the depth of the po- 
tential wells in both are moderately reduced with respect 
to ^^*Hf. This, together with the favorable conditions in 
energy and K value, points to the possibility of very long- 
lived isomers in i86,i88jjf 

The neutron-rich Hf isomeric states are expected to 
decay through f3~ emission competing with 7-ray tran- 
sitions [i. Ill- The K'^ = 8" isomers in i78,i80jjf mainly 



decay to the S"*" member in the ground-state band by 
i^T- forbidden El transitions [l9l |. For neutron- rich nu- 
clei, the 8+ state would lie higher in energy than the iso- 
mer due to reduced /32 deformation. Instead, the isomers 
would decay to the 6'^ member by slower M2/E3 transi- 
tions, making the /3~-decay channel competitive. A sim- 
ilar situation could happen for the 4-qp isomeric states 
where the major 7-transition channel is open to the mem- 
bers in the rotational band built upon the — 8~ iso- 
mer. Spectroscopic measurements may soon be available 
owing to the quick development of experimental tech- 
niques. 

In summary, the multi-qp states in neutron-rich Hf iso- 
topes are investigated by using configuration-constrained 
PES calculations, with attention paid to see the influ- 
ence of deformations on the occurrence of long-lived iso- 
mers. Isomeric states with relatively-low excitation en- 
ergy and very high K value are predicted in Hf isotopes 
at TV w 116. In spite of having small (32 deformations, 
the multi-qp states in -'^^^'^^^Hf remain moderately robust 
against 7 distortion, favoring the formation of long-lived 
isomers. Further work in both experiment and theory is 
needed to study the isomeric structures in neutron-rich 
Hf isotopes, through the N w 116 prolate-oblate shape- 
transition region. 
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